Frictional Instabilities and Carbonation of Basalts Triggered by Injection of Pressurized H2O- and CO2- Rich Fluids

The safe application of geological carbon storage depends also on the seismic hazard associated with fluid injection. In this regard, we performed friction experiments using a rotary shear apparatus on precut basalts with variable degree of hydrothermal alteration by injecting distilled H2O, pure CO2, and H2O + CO2fluid mixtures under temperature, fluid pressure, and stress conditions relevant for large-scale subsurface CO2storage reservoirs. In all experiments, seismic slip was preceded by short-lived slip bursts. Seismic slip occurred at equivalent fluid pressures and normal stresses regardless of the fluid injected and degree of alteration of basalts. Injection of fluids caused also carbonation reactions and crystallization of new dolomite grains in the basalt-hosted faults sheared in H2O + CO2fluid mixtures. Fast mineral carbonation in the experiments might be explained by shear heating during seismic slip, evidencing the high chemical reactivity of basalts to H2O + CO2mixtures

Development and Recovery of Stress-Induced Elastic Anisotropy During Cyclic Loading Experiment on Westerly Granite

International audienceIn the upper crust, where brittle deformation mechanisms dominate, the development of Q4 Q5 crack networks subject to anisotropic stress fields generates stress-induced elastic anisotropy. Here a rock specimen of Westerly granite was submitted to differential stress cycles (i.e., loading and unloading) of increasing amplitudes, up to failure and under upper crustal conditions. Combined records of strains, acoustic emissions, and P and S elastic wave anisotropies demonstrate that increasing differential stress promotes crack opening, sliding, and propagation subparallel to the main compressive stress orientation. However, the significant elastic anisotropies observed during loading (≥20%) almost vanish upon stress removal, demonstrating that in the absence of stress, crack-related elastic anisotropy remains limited (≤10%). As a consequence, (i) crack-related elastic anisotropies measured in the crust will likely be a strong function of the level of differential stress, and consequently (ii) continuous monitoring of elastic wave velocity anisotropy along faults could shed light on the mechanism of stress accumulation during interseismic loading. Plain Language Summary In the upper crust, large strains are accommodated by brittle deformation mechanisms, leading to macroscopic faults embedded within a substantially damaged rock medium. The development of crack damage affects both the strength and the elastic and transport properties of rocks. Nowadays, the evolution of rock elastic properties is commonly used to estimate the direction of the maximum stress along faults and evaluate seismic hazard of seismogenic area. Up to Q6 now, stress-induced anisotropy was expected to be irreversible and observable by geophysics method even after unloading or exhumation of the rocks. In this study, we demonstrate for the first time that unloading induces an almost complete recovery of both stress-induced anisotropy and stress-induced damage. Our results suggest that elastic properties estimated from wave velocity measurement could then underestimate both damage and anisotropy of the crust under shallow depth conditions

Can Precursory Moment Release Scale With Earthquake Magnitude? A View From the Laboratory

Today, earthquake precursors remain debated. While precursory slow slip is an important feature of earthquake nucleation, foreshock sequences are not always observed, and their temporal evolution remains poorly constrained. We report on laboratory earthquakes conducted under upper-crustal stress and fluid pressure conditions. The dynamics of precursors (slip, seismicity, and fault coupling) prior to the mainshock are dramatically affected by slight changes in fault conditions (fluid pressure and slip history). A relationship between precursory moment release and mainshock magnitude is systematically observed, independent of fault conditions. Based on nucleation theory, we derive a semiempirical scaling relationship which explains this trend for laboratory earthquakes. Several natural observations of earthquakes ranging from similar to M-w 6.0-9.0, where precursory moment release could be estimated, seem to follow the extrapolation of the laboratory-derived scaling law. Notwithstanding spatiotemporal complexity in natural seismicity, some moderate to large earthquake magnitudes might be estimated through integrated seismological and geodetic measurements of both seismic and aseismic slips during nucleation. Plain Language Summary Understanding the preparation phase that precedes earthquake ruptures (nucleation) is crucial for seismic hazard assessment because it might provide information on the impending earthquake. Here, we show that the temporal evolution of laboratory earthquake precursors (precursory slow slip, precursory seismicity, and fault coupling) is of little use in assessing an impending earthquake's size. Nevertheless, independent of fault slip behavior (seismic or aseismic) and environmental conditions (stress state and fluid pressure and slip history), the amount of moment released during the preparation phase scales with the earthquake's magnitude. This property is demonstrated by laboratory observations and earthquake nucleation theory and seems compatible with several natural observations of earthquakes ranging from M-w 6.0 to M-w 9.0. As a consequence, if earthquakes exhibit a preparation phase, it could be possible that this phase is larger or longer for higher magnitude earthquakes and consequently, more likely to be detectable

From fault creep to slow and fast earthquakes in carbonates

A major part of the seismicity striking the Mediterranean area and other regions worldwide is hosted in carbonate rocks. Recent examples are the destructive earthquakes of L'Aquila (M-w 6.1) in 2009 and Norcia (M-w 6.5) in 2016 in central Italy. Surprisingly, within this region, fast (approximate to 3 km/s) and destructive seismic ruptures coexist with slow (<= 10 m/s) and nondestructive rupture phenomena. Despite its relevance for seismic hazard studies, the transition from fault creep to slow and fast seismic rupture propagation is still poorly constrained by seismological and laboratory observations. Here, we reproduced in the laboratory the complete spectrum of natural faulting on samples of dolostones representative of the seismogenic layer in the region. The transitions from fault creep to slow ruptures and from slow to fast ruptures were obtained by increasing both confining pressure (P) and temperature (T) up to conditions encountered at 3-5 km depth (i.e., P = 100 MPa and T = 100 degrees C), which corresponds to the hypocentral location of slow earthquake swarms and the onset of seismicity in central Italy. The transition from slow to fast rupture is explained by an increase in the ambient temperature, which enhances the elastic loading stiffness of the fault, i.e., the slip velocities during nucleation, allowing flash weakening and, in turn, the propagation of fast ruptures radiating intense high-frequency seismic waves